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Appearance of Jet-Driving Poynting Flux in Hot, Tenuous Accretion Disks 
Threaded by an Ordered Magnetic Field 
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Abstract 

In a series of our previous works, a model of radiatively inefficient accretion flows (RIAFs) in a global magnetic 
field (so called resistive RJAF model) has proved its ability to account for many physical processes taking place in 
such accretion flows as realized in the nuclei of the galaxies believed to be accreting at a very small fraction of each 
Eddinton accretion rate. Within the present status of this model, however, the model cannot describe the launch of 
a self-confined bipolar jet from the vicinity of disk's inner edge, although it allows the existence of a thermal wind 
widely distributed over the disk surfaces. This is because the electric field (and hence the Poynting flux) vanishes 
everywhere in the disk, whereas such a jet in a globally ordered magnetic field is most likely to be accelerated 
electrodynamically. We show in the present paper that this defect can be overcome naturally if we reformulate the 
problem so as to admit a quasi-stationary change of the magnetic field (and hence the appearance of a non-irrotational 
electric field), and also restore all the terms of order e = (vr/v^) 2 < 1 (where v r and v v denote radial and azimuthal 
components, respectively, of the fluid velocity) which have been neglected altogether in the previous treatments. 
The restored effects are the inertial and magnetic draggings on the infalling matter. As an illustrative example, a 
model solution which is correct up to 0(e) is derived under a set of plausible restrictions. The new solution predicts 
the appearance of a localized Poynting flux in a region near the disk inner edge, suggesting strongly that a jet is 
launched from this region. Another interesting prediction is the appearance of a rapid change of the magnetic field 
also localized to this region. 
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1. Introduction 

Well organized astrophysical jets are often observed around 
various central objects, e.g., both supermassive and stellar 
mass black holes, neutron stars, and even young forming stars. 
Among them, Galactic black-hole X-ray binaries (BHXBs) 
have served as good laboratories for stellar mass black holes 
because of their closeness to the Earth. It has been confirmed 
mainly from their radio and X-ray observations that the pres- 
ence or absence of a jet, and if present, the nature of the jet from 
such a BHXB has close correlation with the state of associated 
accretion disk (e.g., Fender et al. 2004). As often argued by 
many authors (e.g., Fender 2010), such a correlation between 
jets and accretion disks in BHXBs may be considered to apply 
also to their scale-up version, the active galactic nuclei (AGN). 
This conspicuous correlation between the jet state and associ- 
ated accretion disk strongly suggests that the main source of 
the jet-driving power is in the accretion power rather than in 
the spin energy of the central black holes. It is, therefore, natu- 
ral to demand that any accretion disk (or flow) theory that ever 
concerns with the inner regions where a jet might be launched 
should process a due ability to describe this circumstances. 

According to the terminology by Narayan & McClintock 
(2008) inefficiently radiating accretion flows (RIAFs) form 
the low mass-accretion rate branch (see, e.g., Frank et al. 2002; 
Kato et al. 2008) of the series of states called advection domi- 
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nated accretion flows (ADAFs). The other separated branch is 
that for slim disks. RIAFs are specified by a radiation cooling 
time longer than the accretion time, while the slim disks have 
a diffusion time for photons longer than the accretion time. In 
both cases, therefore, the radiation does not have sufficient time 
to carry away significant part of energy from accreting matter, 
and hence main part of energy may be advected down the flow 
in the form of thermal energy. The critical accretion rate which 
divides the two branches seems to be m cr j t ~ 0.01 ~ 0.1, when 
it is normalized to the Eddington rate. 

Currently, the RIAFs are considered to be the most plau- 
sible accretion states that may be accompanied with well de- 
fined jets. Indeed, overall structure of the broad-band spec- 
tral energy distributions (SEDs) from the nuclei of radio-loud 
galaxies, which are the most typical central engines of well- 
developed jets, have been adequately reproduced by the RIAF 
models (e.g., Narayan 2002; Yuan 2007). We know also that 
strong jets appear in the low/hard state of the BHXBs, and the 
state has been understood as the RIAF dominated one (e.g., 
Fender et al. 2004; Narayan 2005). A subclass of the radio- 
loud AGN, the low excitation (or ionization) radio galaxies, 
has been recognized as the counterpart of the low/hard state in 
BHXBs (Baldi et al. 2010) and they also can be understood as 
in RIAF dominated states. 

The basic ideas of RJAF models have been developed by 
many researchers (e.g., Ichimaru 1977; Rees et al. 1982; 
Narayan & Yi 1994; 1995; Abramowicz et al. 1995). One 
of the variations of RIAFs, the convection dominated accre- 
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tion flows (CDAFs), have also been proposed in this context 
(Igumenshchev & Abramowicz 2000; Narayan et al. 2000, 
Quataert & Gruzinov 2000), in which convective transport 
of energy is dominant in the accreting matter instead of adi- 
abatic infall in the case of original RIAFs. A number of re- 
view papers on RIAFs and related flows have been published 
in the literature (e.g., Lasota 1999; Narayan 2002; Narayan & 
McClintock 2008). 

The ability of RIAFs to drive associated outflows has been 
pointed out from the viewpoint of their positive Bernoulli sums 
(i.e., sum of potential energy, kinetic energy and enthalpy) by 
Narayan & Yi (1994; 1995). As a few authors have sug- 
gested (Nakamura 1998; Abramowicz et al. 2000), how- 
ever, the problem seemed to be not so straightforward and 
need much careful considerations. Afterwards, several ana- 
lytical works (Blandford & Begelman 1999; Beckert 2000; 
Turolla & Dullemond 2000) have indicated that what is ac- 
tually needed is to clarify the process of energy transfer from 
the infalling part (i.e., accreting matter) with negative Bernoulli 
sum to the outflowing part (i.e., wind and/or jet) with positive 
sum. Anyway, those outflows modeled explicitly in the above 
mentioned works seem to be widely distributed winds rather 
than much concentrated jets. 

While the gas dynamical version of the RIAF models men- 
tioned above is widely known, it is also widely accepted that 
the presence of an ordered magnetic field in the nuclear region 
is quite essential to the jet formation (Livio 1999; Pelletier 
2004; de Gouveia Dal Pino 2005). Therefore, a model of 
the RIAF penetrated by an ordered magnetic field seems to 
be the most promising version to account for the jet launch- 
ing. Indeed, many numerical works has been performed to 
confirm this point (e.g., Meier et al. 2001; de Gouveia Dal 
Pino 2005). Nevertheless it is strongly desirable to construct 
an analytic model, in order to separately discuss various aspects 
of physical processes taking place in such accretion flows and 
outflows. 

One of such models for RIAF in a global magnetic field 
has been developed by the present author in the last decade 
(Kaburaki 2000; 2001; 2007) and applied successfully to the 
nuclei of some typical dim galaxies (Kino et al. 2000; 
Yamazaki et al. 2002) and radio-loud galaxies (Kaburaki et al. 
2010). Hereafter we call it the 'resistive' RIAF model, in order 
to distinguish it from the well known 'viscous' RIAF model in 
which the presence of a global magnetic field is neglected. 

In the former model, the resistive heating (probably of 
anomalous type) due to the electric current caused in an ac- 
cretion disk plays a dominant role in dissipating the gravita- 
tional energy, instead of the viscous heating in the latter model. 
Although some types of small-scale, turbulent magnetic fields 
are included in the viscous RIAF model through a parameteri- 
zation, it is a great advantage of the resistive RIAF model that 
the model can self-consistently specify the deformation by the 
flow, of an externally given ordered magnetic field. Another 
noticeable difference between the two versions of the models is 
in their disk thickness: resistive RIAFs are geometrically thin 
while viscous RIAFs are thick, although both have a similar 
viral-like high temperature. This is because, in the former case, 
the toroidal magnetic field generated by the action of a rotating 
accretion disk develops mainly outside the disk and compresses 



the disk plasma toward its midplane (Kaburaki 2007). 

Within the resistive RIAF model, it has been shown clearly 
that a thermal wind can be generated according to the state 
of transfer of thermal energy through the plasma (Kaburaki 
2001). The direction (i.e., upward or downward) and strength 
of a wind are specified by a parameter n. When n = the flow 
is adiabatic (i.e., there is no energy transfer through the accret- 
ing plasma) and there is no wind. When n is positive/negative 
(corresponding to the additional input/output of thermal energy 
by some process like conduction, convection or radiation) 2 
there is upward/downward wind from the surfaces of the disk. 
The larger value of \n\ means that the larger fraction of matter 
concerns with a wind. 

Unfortunately, however, the model in the current status can- 
not describe the launch of any electrodynamic jet which is ex- 
pected to occur near the inner edge of an accretion disk. This 
is because the electric field predicted by this model vanishes 
identically in the disk, and hence, also does the Poynting flux 
that is necessary for the jet driving. In the present paper, we 
overcome this defect by improving the accuracy of the model. 

First in section 2, the framework of the resistive RAIF model 
is reconsidered to accommodate the problem of electrodynamic 
jet launching and reformulated as a quasi-stationary problem. 
Then in section 3, each physical quantity is anticipated to be 
written as a multiple of two parts, which are the functions only 
of radius and of co-latitude, respectively, within the geometri- 
cally thin-disk approximation. Such a separation of variables 
is completed in section 4 within a certain degree of approxi- 
mation, and a set of ordinary differential equations is derived 
for functions of radius only. The method adopted in sections 
3 and 4 is essentially the same as in the current version of the 
resistive RIAF model. Section 5 is devoted to the discussion of 
the previously obtained solution, the vanishing Poynting-flux 
(VPF) solution, from our new standpoint. The derivation of a 
model solution, the non-VPF solution, to the newly formulated 
problem is given in section 6. In section 7, the Poynting flux 
and the Bernoulli sum in this solution are calculated. The time 
variability of the magnetic field generated by our model accre- 
tion disks is evaluated also in this section. We summarize the 
results in the final section. 

2. Quasi-Stationary Problem 

As will be shown in section 5, the vanishing of the electric 
field (and hence, of Poynting flux) everywhere in the accretion 
disk in our previous works is closely related to the assumption 
of strict stationarity. Therefore, we need a careful reformula- 
tion of the problem which we really have to deal with. 

For this purpose, we first introduce two characteristic time 
scales associated with the accretion processes in an ordered 
magnetic filed. One is the accretion time defined by 

2 Generally speaking, a non-radiating accretion flow is not necessarily be 
dominated by advection. It may be dominated by convection and/or con- 
duction. Therefore, the name ADAFs should be reserved strictly for the 
special case of general RIAFs (i.e., no-radiating flows), in which there is 
no exchange of energy among the fluid elements (i.e., adiabatic). In that 
case, we can use the term RIAFs for both optically thin and thick branches, 
but have to invent a new name to call optically thin branch distinctively. 
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where L is a representative length scale and v r is the radial 
velocity component (in spherical polar coordinates, r, 6, (p), 
respectively, of the accretion flow. The other is the electrody- 
namic time defined by 

LB 

r cd = — , (2) 

which can be derived from Faraday's law (equation (10) be- 
low). In this paper, we use CGS Gaussian units, and c, E and 
B here denote the light velocity, characteristic sizes of the elec- 
tric and magnetic fields, respectively. 

Quasi-stationary accretion flows in an ordered magnetic field 
are defined as the flows that satisfy the following requirements: 

1. an accretion flow obeys the condition, e = (v r /v v ) 2 < 1, 

2. the fluid motion is stationary, 

3. electrodynamic processes are controlled by the accretion 
processes, and hence, r c d = r acc . 

In the close vicinity of a central black hole, the infall veloc- 
ity may become v r ~ v v or even v r > v v owing to the inward 
attraction of the gravity and effective extraction of angular mo- 
mentum by the magnetic field. Such a flow is usually called 
plunging flow, distinguished from the accretion flows which 
satisfy (w r /w v ) 2 < 1. In this sense, the inner edge of an ac- 
cretion disk (or flow) should be identified with the place where 
v r ~ v v (or e ~ 1). In our current version of the resistive RIAF 
model, the requirement 1 is surely satisfied (Kaburaki 2000, 
2001). 

Note that, in the above definition, the stationarity is required 
only for the fluid motion, but not for the electromagnetic field. 
The time scale for the variation of the electromagnetic field 
is determined automatically from the set of relevant equations 
so as to satisfy the item 3 of the above requirements. In this 
sense, the elctrodynamic processes subordinate the accretion 
processes. Requiring the strict stationarity also for the electric 
field would be too restrictive and constrains the field to be ir- 
rotational (i.e., electrostatic). In this step, the most vital aspect 
of the accretion process may have been lost. 

The requirement 3 above yields 

„ v r B 

E~^—. (3) 

c 

Hereafter, we assume that the accretion flows are non- 
relativistic in the sense S = (v v /c) 2 < («k/c) 2 <C 1, where vk 
denotes the Kepler velocity. This assumption may be justified, 
because the inner boundaries of resistive RIAFs have turned 
out to be located at about a hundred Schwartzschild radii (Kino 
et al. 2000; Yamazaki et al. 2002; Kaburaki et al. 2010), 
where even the Kepler velocity is non-relativistic. Then we 
have 

= £<5<(5«1. 



In a resistive RIAF whose inner edge is located at such radii, 
any general relativistic effect may also be neglected. 

Thus, the set of non-relativistic magnetohydrodynamic 
(MHD) equations applied to the quasi-stationary problem be- 
comes as follows: 



leading equations 
V-B = 0, 

i = f VxB - 

47T 

E=---v xB, 

a c 

v • (H = o, 

(v-V)v = 
P = KpT, 



(4) 
(5) 

(6) 
(7) 



-Vp--^[Bx(VxB)]+g, (8) 
p Airp 



subsidiary equations 

1 dB n „ 
-— = -VxE, 
c at 
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(9) 



(10) 
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The pressure p is assumed to satisfy the ideal gas law, in which 
T and p are the temperature and density of the fluid, respec- 
tively, and K is the gas constant divided by the mean molecu- 
lar weight. The charge and current densities are denoted by q 
and j, respectively. The gravitational force per unit mass, g, is 
simply the Newton gravity which has only the r-component, 

GM 

9r = - — , (12) 

where G is the gravitational constant and M is the mass of a 
central black hole. 

Equation (5) is Ampere's law which neglects the displace- 
ment current in the corresponding Maxell's equation. As stated 
above, the stationarity of the electric field is not automatically 
guaranteed in the quasi-stationary problem. Nevertheless, the 
neglection of the displacement current can be justified since 



1 dE 

c~dt 
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Equation (6) is Ohm's law in which a represents the elec- 
tric conductivity and is assumed to be constant, for simplic- 
ity. Fluid motion is governed by equation of motion (EOM), 
(8), and mass continuity, (7). The time derivative terms in both 
equations have been dropped according to the assumption of 
quasi-stationarity. Also in equation (8), the electric force in the 
Lorentz formula has been neglected since 



\qE\ 







-jxB 

c 







where E/L and j ~ cB/L has been used. Thus, the above 
set of MHD equations are correct in every region of the accre- 
tion flow where e < 1 is satisfied. 

We can evaluate the charge density and the time rate of mag- 
netic filed from the above subsidiary equations, after the elec- 
tromagnetic field has been determined. Local budget of the 
electromagnetic energy is described by the Poynting theorem, 



£~v-*-«-i. 
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which can be derived from Maxwell's equations. Here, u and 
P are the electromagnetic energy density and the Poynting 
flux, respectively, defined by the relations 

u=-L(E 2 +B 2 )*^, P=^ExB. (14) 
«7r 8tt An 

It should be noted here that the energy equation for a fluid 
is lacking in the above set of MHD equations. Usually, this 
equation is replaced by the assumption of polytropic law for 
the fluid, since the inclusion of the former makes the problem 
hopelessly difficult to treat analytically. Also in the resistive 
RIAF model, a parameter analogous to the polytropic index is 
introduced, which is related to the affairs of energy transport 
within it. This point will turn out explicit in section 5. 

3. Axisymmetric, Geometrically Thin Disks 

In the resistive RIAF model, the accretion disk is assumed 
to be axisymmetric and geometrically thin. The former as- 
sumption is only for simplicity but the plausibility of the latter 
assumption is confirmed within the framework of this model. 
Indeed, the model suggests (Kaburaki 2000, 2001, 2007) that 
a strong azimuthal component is generated from a seed mag- 
netic field and develops within the upper and lower halves of 
such a jet bubble as observed in soft X-rays (e.g., Wilson et al. 
2006, Forman et al. 2007, Kraft et al. 2007). The disk is, 
therefore, confined to a thin structure located at the equatorial 
plane of a bipolar jet system, owing to its magnetic pressure. 
When we denote the half opening-angle of an accretion disk by 
A(<C 1), the upper and lower surfaces of the disk is defined as 
6 = ir/2 ± A, respectively. 

As in our previous works, we adopt the method of approx- 
imate separation of variables in order to reach an analytic so- 
lution to the quasi-stationary accretion problem, based on the 
fact that the disk is geometrically thin. The smallness parame- 
ter here is A. Reflecting this fact, it is convenient to introduce 
a normalized colatitude, rj = (8 — vr/2)/A. On the other hand, 
the radius from the center is normalized to the radius of disk's 
inner edge r- m , i.e., £ = r/r in . Note that the definition of £ in 
this paper is different from those in our previous papers: e.g., £ 
in Kaburaki (2000) corresponds to r\ here, and that in Kaburaki 
(2001, 2007) implies £ = r/V out , where r out is the outer edge 
radius. 

Another fact which is useful for simplifying the problem is 
a large deformation of the seed field. In a resistive RIAF, the 
ordered magnetic field penetrating the accretion disk is largely 
amplified from an interstellar seed field by the action of both 
infall and rotational motions in the disk. Since the deformed 
part of the field, b, is much larger than the original seed field, 
Bo, except in the region very close to the outer edge, we use 
the approximation 

B = Ba + b^b. (15) 

Therefore in the set of MHD equations introduced in section 
3, every B can be replaced by b. After this replacement, the 
component expressions of the MHD set in spherical polar co- 
ordinates become to the forms cited in Appendix. 

The seed field B here is considered as a uniform vector 
vertical to an accretion disk, idealizing a large scale interstellar 



magnetic field. The deformed part of the field is maintained by 
the electric current generated in the disk, whose poloidal com- 
ponent has been suggested to close after circulating through 
the boundary of a large scale X-ray cavity and the polar re- 
gions (Kaburaki 2007). The resulting field outside the disk is 
a helical field confined to this X-ray cavity. Such physical con- 
siderations are reflected as boundary conditions to b(£,?y) in se- 
lecting the functional forms given below. Namely, the poloidal 
components of the deformed field b should vanish at large rj, 
while its toroidal component can remain there. The boundary 
condition at the disk outer edge is \b\ — >• |Bo|- We simply re- 
gard that their equality as a definition of the outer edge. 

According to our previous works, we expect separations of 
variables for every relevant physical quantity in the following 



form: 

b r {(,, v) — &r(£) sech 2 ?/ tanh?/, (16) 

bo(t, r))= A6 e (£)sech 2 77, (17) 

r7) = -MOtanh?7, (18) 

jr(£, rj) = -A- 1 ^) sech 2 ?/, (19) 

jete, v)=Je(0 tanhr?, (20) 
3<p(€> r ?) = sech 2 ?? (sech 2 ?7 - 2tanh?y), (21) 

E r te, rj) = A£ r (0 sech 2 77, (22) 

Ee(t v) = ~E g (£,) sech 2 ?? tanhrj, (23) 

E v (£, ??) = -A£^(£) sech 2 ??, (24) 

v r(£, v) = ~v r ({;) (sech 2 ?? — 2tanh 2 7?), (25) 

v$(£, rj) = ASe(£) tanh??, (26) 

*l)=M& (27) 

Pit, V) = P(0 sech 2 ?,, (28) 

V) = P(C sech2? 7> ( 29 ) 

T(£,v7)=T(£). (30) 



Each quantity has been written as a multiplication of its ra- 
dial and angular parts, which are respectively the functions of 
£ and r\ only. The radial-part functions are distinguished from 
the corresponding total functions by tildes on them. The rel- 
ative orders of magnitude among the above quantities are in- 
dicated explicitly in terms of A included in the above expres- 
sions. Therefore, all functions with tilde on them may be con- 
sidered as of order unity with respect to A. It may be inter- 
esting to see that among others v v and T are assumed to be 
independent of 77. 

4. Equations for Radial-Part Functions 

Given the variable separated expressions, (16) through (30), 
the next task is to actually perform the separation of variables 
on each equation cited in Appendix 1 in order to obtain the 
corresponding equation containing only radial-part functions. 
Of course, however, we cannot expect the exact separation of 
variables among so many equations, most of which are so com- 
plex in ry-dependences. Instead, we have to be satisfied with a 
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certain degree of approximation. The status we keep in this pa- 
per (and has been kept in our related previous papers) is as fol- 
lows: i) in equations as many as possible, the separation should 
be completed within the leading order terms in the smallness 
parameter A, and ii) if the separation is not performed even 
among the leading order terms in an equation, it should be done 
at least near the midplane of the accretion disk (i.e., in the limit 
of 77 — > neglecting the terms of order higher than tanh 2 ?y). 

Of course, the separation restricted to the limit of rj — > is 
invalid near and outside the surfaces of an accretion disk (i.e., 
rj 00). In this sense, the above approximation scheme means 
that we focus our primary attention on the main body of the 
disk, where accreting matter is concentrated and expected to 
be governing the physical processes. Another strategy to ana- 
lytically deal with such geometrically thin disks may be the av- 
eraging over the disk height (e.g., Shakura & Sunyaev 1973). 
We believe, however, that the above introduced method of ap- 
proximate separation of variables has many advantages over 
that of height averaging in extracting much physical informa- 
tion contained in the problem. 

Noting that b r ~ 0(1) and bg ~ 0(A), we obtain 



1 9 2 1 db e 



(31) 



as the leading order approximation of equation (Al). The 
term (l/r)dbg/d9, which results from the second derivative 
in equation (Al), is a quantity of order 0(1), since 8/ 89 = 
A~ 1 d/di] ~ 0(A _1 ). The other term (cot 6/r)bg, resulting 
also from the second derivative in equation (Al), has been 
dropped because it is a quantity of 0(A 2 ) in the disk where 
9 — ir/2 < A. Then, the separation can be performed perfectly 
in this leading order equation, yielding a useful expression 



b ^ = L± Hr % r) = i± He ~ br) . 
b r 2dr V ' 2d£, V ^ r ' 

Similarly, equation of mass continuity becomes 



r 2 dr 



{r 2 pv r ) 



1 8_ 

r 89 



(pve 



0. 



(32) 



(33) 



in the leading order approximation, and, further substituted the 
expressions (25), (26) and (29), reduces to 



|£ = r-^-\n(r 2 pv r ) = £— lri(£ 2 pv r ). 
v r dr dt, 



(34) 



Also in this case, the separation of variables is exact among the 
leading order terms. An example of exactly separated case to 
all orders in A is the equation of state, which reduces to 



p = KpT. 



(35) 



The ^-dependences of the current density components, (19) 
~ (21), have been determined from the forms of the magnetic 
field components, (16) ~ (18), through the leading oder ver- 
sions of Ampere's law, 

Att . 1 db 9 

3r = ~ "~ o/i j (36) 
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r 89 ' 
1 8_ 
r dr 
1 db*. 



(rb v ), 



(37) 



This process yields also 

cb 



Anr ' 



38 



Airr dr y v ' 



Cb r 



as a result of variable separations. 

The leading order versions of the component expressions of 
Ohm's law are 

jr 1 

E r = — {veb v -v v bg), 



Eg = --(v v b r 



Jrbip), 



if 1 1 u 

(V r bg 

a c 



vgb r ), 



(40) 
(41) 
(42) 



Both r- and (^-components contain full terms, while the term 
jg j a has been omitted in the ^-component equation. This is 
because we regard A 2 er ~ 0(1) (consistent with the fact that 
plasma is a good conductor), according to our previous works. 
For these equations, the separation of variables cannot be per- 
formed even among the leading order terms so that we must 
be satisfied with the separation in the limit of 77 — > (i.e., 
sech 2 ?/ — > 1, and tanh 2 77 — > 0). The results are 



E r = - u„ 
c 



Jr 

A 2 <7 : 



En 



E a 



(v v b r 



v r bg 



Vrb v ), 



A 2 a' 



(43) 
(44) 

(45) 



Although equation of motion is the most complicated one 
in our set of equations, it can be treated also within the above 
introduced scheme of approximate variable separations. In the 
leading order approximation in A, their three components be- 
come 

,,2 



d_ 
36 



d Vg d 
dr r 89 / 

1 dp GM 

p dr r 2 



ve_d_ 
r 89 



_£_ 

r 



1 



dr 



Anpr 
= 0, 



or 



db r 
'~d9 



(46) 



(47) 



Airpr 



d_ 

dr 



db v 



(48) 



Among them, the ^-component has been so simplified that it 
can be integrated to yield 



v- 



const. 



(49) 



where 'const.' on the right-hand side means that it is inde- 
pendent of 9 (i.e., it may be a function of £ only). When we 
denote this function as p(£), the separation of variables can be 
performed in the limit of 77 — > 0, and we have 

1 ,r, ro. 
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Roughly speaking, this means a balance between the gas pres- 
sure on the midplane of a disk and the magnetic pressure due 
to both <p- and r-components at the disk surfaces. The other 
two equations are separated also in the same limit: 



Kb a 



r 2 



dr Anpr 



1 dp GM 
r p dr r 2 



5 ±(rv = 



(51) 



(52) 



Finally, we go on to the subsidiary equations. The leading 
order approximations of Faraday's law can be written as, 

1 db r _ 1 dE v 
c~8t~ 



r 89 ' 

Idbe Id 

- -Br = - jr\ rE v)' 

c at r or 

1 8\ = _\ (8_ 8E r 
c dt r 1 dr [r e) 89 



(53) 



(54) 



(55) 



The latter two equations are the same as their full versions, i.e., 
no term has been neglected. The former two can be separated 
for any value of rj, resulting in 



1 8b, 
c dt 

1 dbt 



-2^, 
r 

1 d 



c at r dr 
and the latter one results in 



1 dkp 
c dt 



1 



dr 



(rEg 



2E r 



(56) 



(57) 



(58) 



when separated in the limit of 7] —> 0. 

The functional form of the charge density is determined from 
the right-hand side of equation (A16). Keeping only the lead- 
ing order term on the right-hand side, we have 

1 dE e 



H 4nr 89 ' 
and further substituting the expression (23), 

<?(£) v) = — sech 2 ?y (sech 2 ry — 2tanh 2 7y), 



(59) 



Aixr 



(60) 



5. VPF Solution, Revisited 



Under the assumption of strict stationarity, together with a 
few conditions to select a RIAF state, we have obtained an an- 
alytic solution to the set of MHD equations (Kaburaki 2000, 
2001). In that solution, however, the electric field identically 
vanishes in the accretion flow. Therefore, there is no Poynting 
flux to drive such outside systems as bipolar jets often observed 
to be emanating from the nuclear regions. The whole elec- 
tric power generated by plasma motion in the magnetic field is 
merely dissipated locally in the accretion disk, as Ohmic dissi- 
pation. Hereafter, we call this solution the vanishing Poynting- 
flux (VPF) solution. 



In this section, we first observe that the vanishing of the elec- 
tric field is closely related to the assumption of strict stationar- 
ity, and then rederive the VPF solution starting from the van- 
ishment of the electric filed. This section may also serve the 
reader as a brief introduction to the current status of our resis- 
tive RIAF solution. 

It is convenient to rewrite the electric field components, (43) 
- (45), as 



1 



E e =--V r b r (m-^ 



E ip = — v r 
c 



(61) 



(62) 



(63) 



where the poloidal and toroidal magnetic Reynolds numbers 
have been introduced by the relations 



5Rp — 



2 arv r b g _ b v 



(64) 



respectively. The former of the above definitions is introduced 
anew in the present paper. 

The requirement of strict stationarity means the vanishing of 
the left-hand side of equation (10). In our case, this yields 



E v = 0, 



l(rE v ) = 0, 



Er = l±(rEe), (65) 



from its variable separate counterparts, (56), (57) and (58). The 
second of the above equations is always satisfied if the first one 
holds. The implications of the third equation are classified into 
two cases according to the value of Eg. When Eg = 0, the 
equation requires also E r = 0. This is the case of vanishing 
electric field, E = 0. When Eg ^- 0, then dividing both sides 
of the equation by Eg and substituting equation (32), we have 
Eg I rb r = E r /rbg = const.. As will be shown below, the former 
case leads to the VPF solution derived in our previous papers. 
However, we don't know at present whether there is a consis- 
tent set of solution in the latter case. 

Our re-derivation of the VPF solution begins with E v = 0, 
which leads to 5R„ = 1 or 



(66) 



47tA 2 ot L 



The result means that the infall velocity of the accreting matter 
is determined not from EOM but from Ohm's law as a diffusion 
velocity through the ordered poloidal magnetic field. When 



Hp 
result, 



1, both conditions E r = and Eg = yield the same 



(-6- 1 / 2 ) 



(67) 



As stated in section 2, the accretion flow is specified by the 
condition e < 1 so that its inner boundary should be identified 
with the place where e = 1, i.e., 3? 4 = 1 in this case. 

Originally, the VPF solution has been found in the effort 
searching for a state with quasi-Keplerian rotation. The re- 
quirements for a quasi-Keplerian rotation are the following: 
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1. slow infall, i.e., e = (v r /v v ) 2 <C 1, so that equation (51) 
reduces to 



v 2 



1 dp 
p dr 



GM 



(68) 



2. gas temperature with viral form, 

1 dp I GM 
p dr I r 2 



= const., (0<a<l). (69) 



The first requirement is well justified in an accretion disk 
not so close to its inner edge, and when the terms of order e are 
completely neglected, the rotational velocity is determined by 
the gravity partly cancelled by the contribution from the pres- 
sure gradient term. The second requirement assures that the 
temperature is of a virial type (i.e., T oc 1/r) as long as the 
pressure p is a power law function of r. This can be confirmed 
by rewriting it in the form a = KT(—d\np/dr)/(GM/r 2 ). 
From these conditions, we actually have a quasi-Keplerian ro- 
tation, 

i v = (I - a) 1/2 V K , (70) 

where Vk = (GM/r) 1 ' 2 denotes the Kepler velocity. 

The determination of other quantities is straightforward, if 
we introduce the functional form 

b r (0=B in C (3/2 - n \ (71) 

where B- m = Bq £q„ 2 ~™. The form of B; n guarantees that the 
size of the deformed field b should be comparable to that of the 
external field Bq at the outer edge of an accretion disk, £ out . 
The expectation of this kind of power law form for b r comes 
from the fact that T and have already been introduced as 
power law functions. 

The appearance of the unspecified power law index — (3/2 — 
n) reflects the lack of energy equation from our basic equa- 
tions, and the separation of the factor 3/2 is only for conve- 
nience. In this sense, the power law index n may be regarded 
as a free parameter analogous to the polytropic index. Indeed, 
it has been shown that its values reflect the thermal energy bud- 
get of a fluid element in an accretion disk (see, section 7 of 
Kaburaki 2001). Positive and negative values of n correspond 
to a gain and a loss, respectively. The energy transfer may be 
caused by conduction, convection or radiation. Specifically, 
n = corresponds to the case of adiabatic flows, i.e., there is 
no exchange of energy among the neighboring fluid elements. 
Thus, the role of the parameter n is quite analogous to that of 
the polytropic index. 

Substituting the expression (71) into equation (32), we have 
be /K = (2n + l)/4 and hence 

be(t) = ^P-B in C {3/2 ~ n) . (72) 



The result can be used to derive 



«r(0 



r 



(73) 



(2n+l)7rA 2 CTr in ' 

directly from equation (66), and further combined with equa- 
tion (67), 

'•in 



1/2 



Jilt 



(l-a)V 2 V K , 



(74) 



where I4c.in = (GAI/r[ n ) 1/>2 . Simultaneously, we find 5Ri n = 1 
also from the expression (67), since by its definition 5Ri(l) = 1 
should hold at the inner edge (£ = 1) of a disk. Therefore, 
the final expression of the toroidal magnetic Reynolds number 
becomes 

= e 1/2 - 

Then, the result 

MO = B m r (1 - n) 

follows from the definition of 5R t . 

The components of the current density calculated from the 
set of equations (39) are 

m= s ^ L c {2 - n) , (77) 



(75) 



(76) 



cBm 



-(2-n) 



cB h 
47rrj, 



-(5/2-n) 



(78) 



(79) 



Although we can derive the explicit expression for the inner 
edge radius, r in , from the above definitions of v in and $l- m , it 
would be practically useless owing to large ambiguities con- 
tained in the size of the electrical conductivity a. Rather, it 
may be useful to evaluate effective values of the conductivity, 
regarding r- m as a parameter which can be inferred from, e.g., 
spectral fittings (e.g., Kino et al. 2000; Yamazaki et al. 2002). 

Neglecting the term of order e as in equation (68), the gas 
pressure (50) can be approximated like 



Then, substituting equation (76), we have 

p(0 = ^r 2(1 - n) - 

The gas density is calculated from equation (52) to give 

_ 2(n + l) Bf n 



KO = Pin 



<l-2n) 



(80) 



(81) 



(82) 



toril-a) 1 / 2 v in V K ,i 

These results enables us to determine the value of a from the 
definition (69) as 



a = -(1-n), 

where the relations 

1 dp _ 2(l-n)(l-a) 1 /y i y K ,: 
p dr 
GM 



(83) 



^K.in r in 



(2n + l)n, 

r 2 , 



have been used. The allowed range for the parameter n be- 
comes — 1/2 < n < 1 since < a < 1 (in fact, there is a severer 
restriction — 1/4 < n < 1/2, see Kaburaki 2001). 

Substituting the value of a, we obtain the final expressions 
for v r , v v , and p as 

1/2 



Vr(0 



2n+l 



(84) 



>(0 = 



2n + 1 



1/2 



V K - m £ 



-1/2 
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and 



3S? 



;-(l-2t; 



(85) 



(86) 



[Vol. , 
(91) 



respectively. It is also straightforward to derive 



v e (0 = 2n 



and 



no 



l^ 2 



3A' 



2n + l 



r 



1/2 



(88) 



with the aid of equations (34) and (35). The result (87) indi- 
cates clearly that a widely distributed wind (i.e., the vertical 
flow having the same ^-dependence as v r ) appears in propor- 
tional to n, the degree of non-adiabaticity. This is the reason 
why we identify this wind with a thermally driven wind. It is 
interesting to see that also the vertical current is proportional to 
n as seen in equation (78). 

In this section, we have derived the VPF solution starting 
from the condition E = 0, under the restriction of the quasi- 
Keplerian rotation. The derivation is quite straightforward. 
On the other hand, the same solution was derived in Kaburaki 
(2000, 2001) by assuming the condition 



P= i- -r(rv r b v ) I r 2 ^- I — I = const. 
b r dr I ar \ r ' 



(89) 



instead of E r = Eg = 0. Since the condition E v = is a direct 
consequence of db p /dt = (where b p denotes the poloidal 
part of b) and common to both cases, the relation (89) seems 
to be equivalent to the condition E r = Eg = which is the 
one between the two possibilities arising from the condition 
db^/dt = 0. Thus, we have been able to replace the restriction 
(89), whose physical meaning is rather unclear, by the clearcut 
condition E r = Eg = 0. 

Finally, we stress here the advantages of the resistive MHD 
treatments over those of the ideal MHD. In the ideal MHD 
scheme, the electric field is given from the beginning by the 
motional field with negative sign, 



E=--v xB, 



(90) 



neglecting the resistive term in equation (6) as a is sufficiently 
large. However, this asumption poses a very stringent restric- 
tion on the possible solutions. Evidently, we cannot reach in 
the latter scheme the above discussed VPF solution since the 
vanishing of the electric field, E = 0, means v = and/or 
B = as long as the magnetic field is threading the accretion 
disk, which is nonsense. 

In fact, the relative size of the resistive term in Ohm's law 
should be determined self-consistently within a given problem, 
irrespective of the size of the conductivity. Whenever the resis- 
tive term plays a non-negligible contribution, we should inter- 
pret rather that the size of the electric current varies in propor- 
tion to the conductivity, keeping their ratio finite. In the VPF 
solution, we have 



v x B, 

a c 

since E = 0. This means a local balance of the motional field 
with the resistive term. Therefore, the electromotive force does 
not remain outside the driving system, the accretion disk. The 
whole generated electromagnetic energy is dissipated locally 
as the Joule heating. We stress again that this kind of solution 
can never been obtained in the ideal MHD scheme. 



(87) 6. Non-VPF Solution 



In the previous section, we have seen that requiring a station- 
ary magnetic field is closely related to the vanishment of the 
electric field, at least, under the assumption of quasi-Keplerian 
rotation. Since such a solution cannot explain the generation 
of jet-driving Poynting flux in an accretion disk, we have to 
abandon this assumption in order to proceed toward jet launch- 
ing. For this reason, we have formulated the quasi-stationary 
problem in section 2. At the same time, we have to deal with 
the inner region of the accretion disks more carefully because 
jets are known to be launched from the vicinity of disk's inner 
edges, where e approaches unity. This means that the complete 
neglection of 0(e) terms is insufficient. 

Thus, we now return to the original set of MHD equations 
obtained in section 4, without neglecting the left-hand side of 
equation (68) and the magnetic pressure term arising from b r in 
equation (50). In this new framework, we try to extend the VPF 
solution to be correct up to the first order in e. Of course, even 
the first order corrections to the VPF solution are not sufficient 
to describe exactly the regions very close to the inner edge and 
below it. Nevertheless, it is well expected that the result may 
give us a certain clue to the physical origin of jet-driving pow- 
ers. 

In order to find a model case of such a solution (hereafter we 
call it simply the non-VPF solution), we put a few restrictions 
on the solution: 

1 . the non-VPF solution should coincide with the VPF so- 
lution in the limit of £ — > oo (i.e., e — > 0), 

2. as an extension of equation (69), we require 




1 



n).(92) 



The requirement 1 is quite natural, and serves as a boundary 
condition for each physical quantity. In accordance with this 
requirement, we fix hereafter the position of disk's inner edge 
by the relation eo = 1, where eq is the value of e calculated from 
the VPF solution, i.e., eo = ("r/^JvpF = Therefore, this 
place can be different from that calculated from e = 1 in the 
non-VPF solution. 

On the other hand, the requirement 2 may seem somewhat 
intuitive or arbitrary. Although it is easy to confirm that this 
condition is identical with (69) in the limit of £ — > oo, one may 
feel that there is no necessity to require (92) among other such 
possibilities. We agree with this assertion, but it is also true that 
if we admit this condition we can reach a set of physical quan- 
tities all of which satisfy every equation within the required 
accuracy. Therefore, it can be called a solution. Of course, 
however, we never insist that this is the unique solution. In this 
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sense, the non-VPF solution obtained in this paper should be 
called a model solution. 

Eliminating the pressure p in equation (51) with the aid of 
equation (50), we can rewrite the r-component of EOM as 

_ dv r bl 



K be d 

— '— r—\nb 

dr Airpr \ b r dr 



v% 1 d 



r p dr \ 8ir 




GM 



(93) 



where the terms of 0(e) have been collected on the left-hand 
side. Further substituting equations (32) and (92), we finally 
obtain 



r— lnu r 
dr 



2D dr 




(94) 



where the discriminator of radially infalling matter have been 
defined by 



D 



(95) 



with Va = 6 r /(47T/5) 1//2 being the Alfven velocity. 

Since the factor v 2 r on the left-hand side of equation (94) 
is already a quantity of 0(e) by definition, it is sufficient to 
evaluate the factor in the curly brackets by using the lowest 
order solution, i.e., the VPF solution. The values of the related 
quantities evaluated in terms of the VPF solution are 

[£]vpf = ^^, (96) 



d r d ( b r 

r d-r lnVr+ 2Dd-r ln [^ 



2?i + 9 
' 2n+T ' 



(97) 



J VPF 

It is evident from the latter result that the left-hand side of 
equation (94) is negative, and hence, it follows that v v < 
(1 -aY^VK. This means that, owing to the inclusion of the 
inertial and magnetic draggings on the left-hand side, the cen- 
trifugal force required to balance the gravitational attraction 
can be reduced from the VPF case (especially in the inner re- 
gions). In the limit of > oo, we have v v — > (1 — a) 1 / 2 Vk from 
equation (94), neglecting the terms of ~ (e) on the left-hand 
side. 

Although these constraints are not adequate to determine the 
functional form of v v uniquely, here we anticipate its from as 

1/2 



,(0 



2n + 1 



V K ,in £ 



-1/2 -AC 



(98) 



in order to obtain an illustrative model solution of the non-VPF 
version. Substituting this expression into equation (94), we 
obtain the expression for S 2 containing a parameter A. The 
value of A can be determined by the condition, v r — > [u r ]vPF 
as £ — > oo. The results are 



5r(0 



2n + l 



1/2 



v KM c 1/2 



i 



-2AC 



1/2 



2.4 



2n + l 



1/2 



^inf 1 (for£»l), 



where 



A = 



2n + 9 
12 



(99) 



(100) 



The value 1 — a = (2n + 1) /3 has been substituted in the above 
derivation. 

The last expression for u r (£) is obtained by expanding the 
exponential factor in a power series of eo = and keep- 
ing only the terms up to the first order. In the accretion disk 
(where £ > 1), the approximation to this order is sufficient for 
the present purpose. Namely, the inclusion of the first order 
corrections in e has not altered the form of v r from that in the 
VPF solution, equation (84). It should be noted, however, that 
the former has been derived from the r-component of EOM in- 
stead of Ohm's law. The ratio of the velocity components then 
becomes 

-1/2 



1 



-2A£~ 



2A 



1/2 -AfT 1 



(101) 



This means that e = (vr/v v ) 2 ps exp(2A£ _1 ). Thus, the 
actual value of e deviates largely from eo near and below the 
inner edge (£ = 1). 

Next, we rewrite the ^-component of EOM, (52), with the 
aid of equation (95). Eliminating 1/47T/5 from the former, we 
obtain first 



bebif 
bl 



D 



dl/dr 



1 = 



and further using equation (32) and the definition of 5R t 



dr 



r— ln(r 2 6 r ^ 
dr 



(102) 



Similarly, the extended RIAF condition (92) can be rewritten 
in the from 

2 



-r— ln& 

dr 



2(1 -n) 



<p 



D 



(103) 



2n + l \Jft t v ri 

Since the ratio v v /v r in our model case is approximated by the 
last expression in (101), the above equation reads 

d ~ 1 



-r— \n.b v 
dr 



(104) 



As a trial function for b r , we chose here 

b r (0 = Bin £- (3/2-«) e -(2n+DA«- 1 . (105 ) 

The exponential factor added to the corresponding VPF solu- 
tion causes a rapid decrease of this component at radii smaller 
than the disk inner edge. This seems very reasonable since it 
should vanish on the polar axes owing to the symmetry require- 
ment. The above b r results in 



MO 



2n+l 



X (1 + 2AC 1 ) £-(3/2-n) e -(2«+l)A ? - 1 ; (1Q6) 

from the flux conservation equation (32). Then, further noting 
that rd\nl/dr = (1 + 2A£~ 1 )/2, we obtain 
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2n+l 

from equation (102). 

It is very interesting to examine the possibility in which 

2n + l 



(107) 



D(0 



(108) 



Of course, this satisfies the requirement D — > (2n + 1) /2 since 



— > 1 in the limit of £ — > oo. The discriminator selected 



above corresponds to a trans-critically accelerated infall for any 
allowed value of n (— 1/4 < n < 1/2). Such an infall is analo- 
gous to that in Bondi's theory of spherical accretion, although 
no special affairs appear at the 'critical point' (i.e., at the point 
where v r = Va or D = 1) in our problem. If the discriminant 
is given by the above formula, then equations (107) and (104) 
reduce, respectively, to 



and 



-r— In b. 

dr 



(1-n) e~ M \ 



(109) 



(HO) 



The form of 5Rt is thus identical to that in the VPF solution. 
We can easily confirm that the functional form 



-(X-rO.-Cl-nJAr 1 



(111) 



satisfies equation (110) within the required accuracy, since 1 — 
j4£ _1 w e~ A ^ up to £ _1 (= e o)- In order for this b v to be a 
non-VPF solution, it should be consistent with equation (109). 
Since we have 

M t = ^=e /2 e 3nAr \ 

b r 

it turns out that the three components of the magnetic field in- 
troduced in this section and the postulate (108) are consistent 
with the set of quasi-stationary MHD equations if and only if 
n = 0, i.e., the accretion flow is adiabatic. Thus, the value of 
the parameter n, which cannot be determined in the framework 
of the lowest order approximation, is restricted to the special 
case of n = in our improved treatment, without referring to 
the energy equation. This fact suggests that the realization of 
the present, non-VPR model solution in actual situations sensi- 
tively depends on the affairs of energy transport which can be 
realized in the accretion flows. 

Thus, the final forms of the magnetic field components in 
our model solution (in which n = 0) become 

M0=A n r 3/2 e- Ar \ 



1 



-l V -3/2 e -Ar 1 



b e (0 = -B in (l + 2AC L )Z 

MO = -Bin r'e-^ 1 , 

with 

The discriminator of the radial flow is specified as 
1 



D(0 



(112) 
(113) 
(114) 

(115) 
(116) 



The derivation of remaining quantities are straightforward. 
We obtain 



SB? 



from equation (52), 
from equation (50), and 

r(0 = ^(i+r 1 )' 



(117) 



(118) 



(119) 



from equation (35). Here, the deep suppression of the temper- 
ature from the corresponding VPF value near the inner edge 
should be noted, since it has already been expected from spec- 
tral fittings for some LLAGN based on the VPF solution (e.g., 
Yamazaki et al. 2002; Kaburaki et al. 2010). The suppressions 
in the density and pressure are consistent with the presence of 
outflows from the accretion disk in the forms of jets. 
The final (i.e., when n = 0) forms of v r and v v are 

5r (0 = %r 1 , (120) 



V3 



V3 

and equation (34) yields 
AV K 



M0 



V3 



(121) 



(122) 



As confirmed from the angular dependence in (26), v$(£, rj) 
resulting from the above vg (> 0) describes an outflow from 
the accretion disk. Different from the case of the VPF solution, 
vg here exists even in the case of n = 0. Corresponding to this 
fact, the right-hand side of equation (122) contains the constant 
A instead of n (compare with vg(^) in the VPF solution). The 
functional form of this outflow is also a power law but steeper 
than the radial infall of the present solution (see figure 1) and 
the vertical outflow in the VPF solution. These are the reasons 
why we expect that the vertical outflow in the non-VPF solu- 
tion corresponds to the base of an electromagnetically driven 
jet, distinguished from the thermally driven wind in the VPF 
solution. 

7. Poynting Flux, Time Variability and Bernoulli Sum 

The components of the current density in our non-VPF 
model solution are derived from Ampere's law as 

~m = ^r 2 e- A ^\ (123) 



MO 



U(0 



cAB u 
4.7rnu 
cB hl 
47rr in 



-5/2 p -AV 



(124) 



(125) 



Similarly to the case of vg, jg contains the constant A instead 
of n in the corresponding VPF versions. The presence of this 
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Fig. 1. Radial profiles of the 9-component of the Poynting flux (Pg) 
and the velocity components of the accretion flow (v r , vq, v v ). They 
are normalized to AVK,mB? /4V37T and VK,in/ v3> respectively. 
The profiles are written based on the functional forms obtained in the 
text and shown extended to the inner region (£ < 1) where the accu- 
racy of the solution may be poor (v r is based on the first expression 
in equation (99)). The downward deviation of the rotational velocity 
{v v ) from the Kepler curve becomes appreciable below £ ~ 3 and the 
vertical velocity exceeds the radial one only below £ ~ 1, The vertical 
Poynting flux has a sharp peak near the disk inner edge (eo = £ = 1). 

component even in the case of n = would be important from 
the viewpoint of current closure. 

In order to calculate the full components of the electric field 
from Ohm's law, we need the value of 3? p . This can be evalu- 
ated from the definitions of V[ n in (73) and of 5R P in (64) as 



®p(0 



1 + 2AC 1 , 



(126) 



271+1 Wi n b r 

where the value of u; n has been fixed by the condition = 1. 
Using this result, we obtain 



Er(0 - 



4\/3c 

VK,inAr 



V3c 

A Vk, in B in 

2\/3c 



(127) 



1-e 



-Ar 1 



(128) 



(129) 



The calculation of the Poynting flux is straightforward, and 
the resulting components are 

(130) 



p (n _ AV KM B in 4 2AZ- 1 



16V^7T 
^4VK,inS; 2 



(i+2r 



-9/2.-2^" 



3/2^ (131) 



(132) 



4^7T 

The existence of a steep negative power-law dependence in £ 
and that of an exponential factor in each of these components 



clearly indicates a localized nature of the Poynting flux to the 
region near or even within the disk inner edge (see figure 1). 
Since we have 



V-P = 



P0® 



sech 77, 



(133) 



the positivity of P§ guarantees that the accretion disk is an elec- 
tric power source and is driving outside system, most probably, 
a bipolar-jet system. Although the structures of the resistive 
MHD jets at rather remote regions from their footpoints have 
been discussed in our previous work (Kaburaki 2009), fur- 
ther detailed investigations are required to fully describe the jet 
launching processes. 

In our quasi- stationary treatment, the time variabilities in the 
magnetic field can be calculated according to equations (56) ~ 
(58) as 



ldb r 
Vdt 

c dt 



AVkahBu 
2^/3 cr 
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4^ 



cr 



-9/2 e -A£~ 



i-9/2 £ -AC 



(134) 
,(135) 



1 8b v AV K ,inB in 



dt 



2y/3 cr 

l_ 4r V2 



2 AS, 



-3/2 
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(136) 



As is evident from the presences of the exponential and steep 
power-law factors also in these expressions, the time variations 
become noticeable only in a narrow region near the disk inner 
edge. The characteristic time scale of these changes is that 
of the infall according to the requirement 3 in section 2, and 
is very short in this inner region supposed to be the base of 
a jet. This fact is very interesting because such a rapid time 
variation in the magnetic field may well explain the presence of 
rapid phenomina such as flares or bursts often observed at the 
roots of various jets. For example, we can imagine that a large 
accumulation of the magnetic field would cause an impulsive 
acceleration of fast jet. 

The Bernoulli sum is the sum of energy carried along with a 
fluid element per unit mass, and is defined here as 



1 



B e (0 = -(v 2 r + v 2 J 



5 p 
2 ~P 



GM 



(137) 



where the second term on the right-hand side is the enthalpy of 
an ideal gas. The kinetic energy associated with vg has been 
neglected because it is a small quantity of C(A 2 ). Here, we 
pay attention only to B c (£) although the 77-dependence of each 
term is different from each other. This is because our primary 
interest is on the midplane of an accretion disk around which 
the majority of matter is concentrated. The lack of the electro- 
magnetic field energy in the sum is due to the fact that the field 
is distributed in space and is not carried by the fluid motion 
in the quasi-stationary situation (in contrast to the ideal MHD 
case). 

For a discussion of overall behaviors of the Bernoulli sum, 
we have to extrapolate the results of our non-VPF model so- 
lution beyond the inner edge toward the center (i.e., £ < 1). 
For this purpose, the original expression for iy given in (99) 
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is suited. The final result calculated in terms of the non-VPF 
solution is 



- v. 



K. 



i- i 



15 

Ic 



-(3/2)r 



(138) 



where the value A = 3/4 has been substituted. The asymptotic 
values for large and small £ are 

1 



9 K K,iB? , 



as £ —> oo, 
as £^0. 



(139) 



In the VPF solution, we have B e (£) = within the accretion 
disk (£ > 1) when n = (Kaburaki 2000; 2001). The above ex- 
pression for the limit of £ — > oo includes the first order correc- 
tion in to that result. In both large and small radius limits, 
_B (£) behaves like £~ 2 having negative and increasingly large 
absolute values toward the center, which reflects the nature of 
infalling matter. 

8. Summary and Conclusion 

The resistive RJAF (radiatively inefficient accretion flow) 
model has been very successful as a model for the inefficiently 
radiating accretion disks threaded by an external, ordered mag- 
netic field. It provides a firm basis to discuss various physical 
processes taking place in such accretion flows. It can also well 
reproduce the main properties of observed broad-band spectral 
energy distributions emanating from the nuclei of nearby nor- 
mal galaxies and some types of active galactic nuclei, which 
are all believed to be accreting at a very small fraction of each 
Eddington accretion rate. The present status of this model, 
however, is unsatisfactory in that it cannot explain the electro- 
dynamic formation and acceleration of bipolar jets which are 
frequently observed associated with such accretion flows. This 
is because the Poynting flux which is necessary to drive outside 
systems (especially a jet) vanishes identically. 

In order to overcome this defect in the current version of our 
model, we have carefully reformulated in the present paper the 
problem we have really to deal with. The new scheme allows 
a quasi-stationary change of the magnetic field, and hence, the 
presence of a non-irrotational electric field. Another imrpove- 
ment in the new scheme is the restoration of the terms of 0(e), 
where e = (v r /v v ) 2 , which have been neglected as small. This 
is to keep the problem valid even in the inner regions where 
e ~ 1. The main restored terms are the inertial force of the in- 
falling matter and the dragging force due to the poloidal mag- 
netic field. Both these terms appear in the radial component of 
EOM, and reduces the rotational velocity of an accretion disk 
near its inner edge. 

As an example of the solutions to this modified problem, the 
non-VPF (vanishing Poynting-flux) solution which is accurate 
up tp the first order in e, has been derived under some plausible 
constraints. Although this solution has many essential advan- 
tages over the previous version (the VPF solution), it is still an 
approximate solution in which the accuracy becomes poor near 
the disk inner edge. Therefore, one should be well aware of the 



possible inaccuracies when he applies the obtained results to 
such regions. 

The obtained solution has a non-vanishing electric field and 
a diverging Poynting flux, the latter of which suggests that the 
accretion disk is electrodynamically driving some outside sys- 
tems. Interestingly, the model holds only when the RIAF is, in 
particular, an ADAF (advection dominated accretion flow), that 
is, when the accretion flow is adiabatic and hence no thermally 
driven wind from the disk surfaces exists. The predicted verti- 
cal components of the Poynting flux and the outflow at the disk 
surface are much localized to the region very close to the disk 
inner edge, strongly suggesting that the resistive RIAF disk has 
an ability to drive electromagmetic jet. This result seems to be 
quite consistent with the observational fact that strong jets are 
often associated with the nuclei accreting at very small frac- 
tions of their Eddington rates. In such rarefied accretion disks, 
the radiation losses and energy exchange among accreting mat- 
ter cannot occur efficiently. 

The obtained model also suggests the presence of rapid 
changes in the magnetic field localized to the inner accelera- 
tion region. This fact may explain the origin of some violent 
phenomena such as the flares and/or impulsive acceleration of 
strong jets. 

Thus, our resistive RIAF (or ADAF) model has been greatly 
improved to be able to admit the formation of electrodynamic 
jets in the region very close to the inner edge of an accre- 
tion disk. However, for the detailed understanding of con- 
crete mechanisms of the jet launching further investigations are 
needed. 

Appendix 1. Expressions in Spherical Polar Coordinates 

The expressions in spherical polar coordinates of the basic 
equations for the quasi-stationary problems are summarized 
below. The assumption of axisymmetry has already been taken 
into account. 



• magnetic flux conservation: 
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• mass continuity: 
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• equation of motion: 



(smdpv e ) = 0, (A8) 



d vg d 
dr r 88 / 
1 dp GAI 
p dr r 2 

d_ 

dr 1 



1 



i 



Airpr 
db, 
~d9 



, (A9) 



d v e d \ 1 , 2 

«r^-H TTTi Ufl + -(fr«e-'" ¥ ,COt< 

dr r dO I r v 



1 dp 1 
pr 90 47rpr 



Sin 6/ W 



, (A10) 



dr 



— "ht; hv-i (« r + wecot( 

r d0 I r 



1 



4-7rpr 
equation of state: 
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Faraday's law: 
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• charge density: 
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